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Abstract 

We complement our earlier position-space exploration of a recently proposed S-matrix 
for transplanckian scattering by a momentum-space analysis. As in the previous paper, 
we restrict ourselves to the case of axisymmetric collisions of extended sources. Com- 
parison between the two formulations allows for several cross-checks while showing their 
complementary advantages. In particular, the momentum-space formulation leads to an 
easier computation of the emitted-graviton spectra and to an attempt to study the system 
beyond its critical points into the presumed gravitational- collapse regime. 
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1 Introduction 



In a companion paper [T] we have explored a recent proposal [2j for an S'-matrix descrip- 
tion of transplanckian scattering in four-dimensional spacetime. The proposal is based 
on an approximate resummation of the semi-classical corrections to the leading eikonal 
approximation [3] and amounts to solving the equations of motion of an effective action 
introduced quite sometime ago [HIS]. In a suitable limit, the longitudinal dynamics can 
be factored out leaving behind an effective dynamics in the "transverse" two-dimensional 
space. 

In [1] we considered, quite systematically, the case of the axisymmetric (central) col- 
lision of two extended sources/beams. This case has several advantaged over the original 
problem of a two-particle collision at non-vanishing impact parameter: 

• The partial differential equations (PDE's) reduce to ordinary differential equations 
(ODE's) making the problem affordable by analytic and numerical techniques, with- 
out having to make either an azimuthal-average approximation [2] or to resort to 
advanced numerical techniques [B] . 

• The IR-sensitive graviton polarization, which was neglected by hand in [2j , is simply 
not produced in the axisymmetic case [1]. 

• We can consider a wide variety of initial states by playing with the many (shape 
and intensity) parameters chracterizing the sources and check for the existence of 
critical surfaces in this multidimensional space. The results can then be compared 
with those coming from closed trapped surface (CTS) criteria [7], [8] and will be 
tested, hopefully in the near future, against numerical GR calculations (see [H] , [ID] 
for a few results already available for this case). 

The results of [1], based on a position-space analysis, gave further support to the con- 
clusions of [2]. In particular, we were able to prove a one-way relation between the CTS 
criterion of [5] and the criticality condition in the ODE system. We could also determine 
quite precisely the critical surfaces in a variety of cases and found a good quantitative 
agreement with CTS-based criteria [8]. Finally, we confirmed that, above those criti- 
cal lines, a new absorption of the elastic S-matrix turns on with a universal behaviour 
reminiscent of Choptuik's critical exponent pT] . 

In this paper we complement our previous work [T] by a momentum-space analysis 
according to the following outline: in section 2 we present the momentum space formula- 
tion of the extended-source problem and, after specializing to the axisymmetric case, we 
give the explicit form of the action and of the equations of motion. In section 3 we recall 
some interesting extended sources already considered in [1] adding their momentum-space 
form. In section 4 we study numerically the field equations, determine the critical lines 
in parameter space, and compare them with the position-space results of [T]. In section 
5 we discuss the spectrum of the emitted gravitons starting from the perturbative regime 
and until one approaches the critical lines. In section 6 we present an attempt to extend 
the solutions in the (presumed) BH-phase. 
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2 Momentum space action and field equations 

We recall from [Ij the position-space action of [2j generalized to extended sources: 



A 



-ViaVia 



a[x)s[x) + a[x)s{x) 
£x ((VV)^ + 20(V^a V^a - ViV^a ViVj-a)) 



with the three real fields a, a and (j) representing the two longitudinal and the single, 
IR-safe, transverse component of the gravitational field, respectively. 

The center of mass energy y/s provides the overall normalization factor 2txGs = ^R"^, 
while the two sources s{x), s{x) are normalized by / d'^x s{x) = J d'^x s{x) = 1. In order 
to go to momentum space we can either start from ([1]) or generalize directly eq. (5.2) of 
[2]. The result is: 
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-h{k)h{-k) - h{-k)n{k) 
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where the FT of the sources (still denoted by Sj) are normalized by requiring Sj(0) = 1. 
Furthermore, 



n{k)^ 



, p,{k) = ^^ , h{k) = -k'm 



d'^kid'^k25{k -ki- fc2)/9i(fci)/32(fc2)W^i2 



^2^2 



(3) 
(4) 



and 9i2 is the angle between the two transverse momenta ki and ^2. 

If we now specialize to the axisymmetric case where sources and fields depend only 
upon kf, we can use the following relations: 



d k = 7T dk 



d'^kid'^k25{k — ki — k^ 
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^ J\>0 l^l'^2'^ 



s^n^e^2 1 f dkjdkj /;^(^2^^2^p) 



where 



and we have used the identity: 
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together with the change of variable (at fixed k"^, kf) 

dkl = d{k — kiY = 2|A;| |A;i| sin OkidOki = 2\k\i\k2\ sin OudOki ■ (8) 
One thus arrives at the final form of the axially-symmetric effective action: 

^ = J^[m^s2ik') + /3,ie)s^ie)-Me)^,ie)] 
^""^^^ jdk^\h{k')h{k') + ^ j '^^l^^f' ^xiki ki k^)h{e)mimkl) , (9) 
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whose equations of motion read 
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A(fc^) = s,{k') + -^\lKkl kl k^)h{klmkl) ■ (10) 

It is not a completely trivial exercise to show directly the equivalence of these equations 
with the corresponding ones in position-space p| (where a dot stands for d/dr"^), 

2Tip{r) R ' 
l.o_.^2. lRi{r)R2{r) 



p = -{2'kR) (21(22 = - 



2' ' 2 p\r) ' 

i?.(r) = R I d^x Si{x) . (11) 

itx|2<r2 

The proof, not reported here, makes use of the following (known?) integral of three Bessel 
functions (that we have checked numerically): 



' drJ,{rk)J,{rk,)Urk2) = ^ VKKkuk,) ^^^^ ^ ^^2) 

Zn KK1K2 

3 Examples of source profiles 

In this section we list various extended sources already introduced in [1] and give their 
momentum representations. 

A. As a first class, consider finite-size sources with the following profiles: 

Later, without lack of generality, we shall be fixing the transverse size of the two identical 
beams L to be 1. One can easily verify that these sources satisfy our normalizations and 
that 



/ dph{p) = R{r/L)/R=—- ^ s{r)dr^ = l. (14) 

Jo {L^d + r\l - d)) ' Jo 
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Once Fourier transformed to momentum space and normalized according to the pre- 
scription of section 2, the above sources become {k = |k|): 

Si(fc) = S2{k) = p2 = / dr^Jo{kr) -r- . (15) 

In particular, for two homogeneous beams {d = 1) we have: 

si{k) = S2{k) = j dr'Mkr) = -^MkL) . (16) 

B. Point-like sources are difficult to deal with numerically, especially in momentum 
space. We introduce therefore Gaussian-smeared versions of the point and ring-like sources 
considered in P: 

.i(x) = exp (-^) ( ^1 - , s^ix) = exp (-^^^^^) ©(^2 - r). 



M ^ - exp (gl) ) , AT, ^ 2. (.^(exp (0) - 1) + ^^./|Erf-^) 



(17) 



When 0" — > (oo) such configuration reduces to the one of the point-ring (two homoge- 
neous beams) case. The corresponding Fourier transforms are: 

Si{k) = 271 I rdrJo{kr)si{x) . (18) 
Jo 

C. Another interesting example is that of gaussian sources concentrated at r = 0. 
They correspond to: 

or, in momentum space, to: 

/ u-2 r 2 \ 

s,ik) = exp . (20) 



4 Numerical solutions and comparison with x-space 
results 

There are two ways to solve the non-linear system (fTOll . One may use an iterative(recursive) 
procedure suggested by the form of the equations, or treat them (after discretization) as 
an algebraic system of polynomial equations of third order. The two approaches are to 
some extent complementary. The recursion turns out to be convergent only in the disper- 
sive phase, and can therefore be used to determine the inter-phase boundary in parameter 
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space. The algebraic approach allows to explore also the BH phase by generating genuine 
complex solutions of the system. It can also be cross checked, of course, with the recursive 
method in the dispersive phase. 

Both approaches rely on a momentum discretization procedure. As a first step, we 
introduce new variables which span unit intervals and are thus convenient for that purpose: 

""'i^hL' y~YThL' """T+fcL- ^^^^ 

Here L is the size of the two identical sources (A, B from the previous Section) or of one 
of them (Li in the case C). In the following we set L = 1. In the new variables equations 
(fTOl) read: 

TT^ x{1 - x)y{l - (1 - vy 
Piiv) = s^iv) + ^ ^ 3.(1 J^^yli _ v)h{x)/3,iy) . (22) 

At fixed V the x, y integrals are over the triangular region T which is bounded by three 
hyperbolas 

0<t;<l, 0<a;<l, y.minix, v) < y < ymax{x, v), (23) 

f \ ( \ i V-T+XV ^ < fr)A\ 

X + V-XV { X>V 

Next, we discretize the variables 

1 i — I 

u — >Ui = —^ , i = l, ...,n, u = x,y,v, (25) 

2n n 

and turn the integral equations into a set of 3n algebraic equations 



K 
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= + :-S;;,,=i^g J2n+./. , (26) 
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2 



fn+i — Sn+i + — S"^fc=iW-J^^/2n+i/n+i , (27) 
f2n+i = '^^],k=l'^i,j,kfifn+i ; (28) 

where the weights w are the discretized versions of the kernels * measure in the corre- 
sponding continuous equations 

(31) 
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and the Area{i, j, k) is the area of the intersection of a small, 1/n^, square with the 
"triangle" T if the center of the square {xj,yk) is inside T, and zero if the center is 
outside T. 

As already remarked, these equations can be used either to set up a recursive proce- 
dure, or can be directly solved numerically. The latter approach, dubbed the "algebraic 
method" , works in both the BH and the dispersion phase. 

In the first two (A and B) rows of Table [T] we compare a sample of critical values of 
R (in units of L) as determined by configuration [1] and momentum space methods. In 
the latter Rc is defined as a point where the recursion (126!) diverges. As such it depends 
on the length of a trajectory n and in principle requires extrapolation to n = oo. In the 
Table we used n < 80. 



d 




0.5 


1.0 


1.6 


2.5 


4.0 




A-x 




0.419 


0.471 


0.502 


0.528 


0.550 




A-p 




0.429 


0.476 


0.499 


0.501 


0.477 




a 




0.01 


0.1 


0.2 


0.3 


3.0 




B-x 




0.615 


0.572 


0.525 


0.486 


0.470 




B-p 




0.058 


0.436 


0.501 


0.489 


0.476 




P 


0.25 


0.333 


0.5 


1.0 


2.0 


3.0 


4.0 


C-x 


.810 


.816 


.821 


.823 


.821 


.816 


.810 


C-p 


.823 


.833 


.850 


.841 


.838 


.840 


.832 



Table 1: {R/L)c for a range of sizes of the power-like and Gaussian sources: a comparison 
between configuration and momentum-space results. A, B and C label sources as discussed 
in Sect. 3. In the case C: p = L2/L1 and the critical value of the ratio 2R/{Li + L2) is 
shown. 

For extended sources both approaches are consistent meaning that at n ~ 50 — 80 
momentum space estimates have already converged. However, for narrower sources, mo- 
mentum method requires a yet finer discretization. This is to be expected, since a finite 
mesh in momentum, say Aj9, limits the spatial resolution to Ax > 1/Ap. In such cases 
one has to extrapolate numerical data from the case of extended sources as done in [6]. 
In the following, we shall be discussing only homogenous beams where a finite n ~ 50 is 
adequate. Notice again the two special cases mentioned earlier, namely d = l(oo) and 
(J = 00(0) which correspond to the scattering of homogeneous beams and that of a par- 
ticle and a ring. The critical radii for these cases (Rc ~ .47 and Rc = 2^^'^3^^^^ ~ .62 
respectively) were obtained in [2.J and agree with the ones quoted in the Table except for 
the momentum study of the infinitely narrow sources which was to be expected. 

The third row of the Table summarizes the head-on collision of the two central, gaus- 
sian sources with different widths (C). The problem is symmetric with respect to the 
interchange Li ^ L2 therefore we display Rc in units of (Li +L2)/2. Agreement between 
X- and p-space methods is quite satisfactory. The configuration space technique used in 
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pr| was manifestly symmetric under the interchange of Li and L2 as reflected in the Table. 
However, in the momentum space calculations we have deliberately used only one source 
size as a scale. The resulting small asymmetry gives an idea of the sensitivity to the 
discretization parameter n (which was not so large here, n < 20). 

Finally, we emphasize a very weak dependence of [2R/{Li + L2)]c on p. This confirms 
the observation, made already in pQ, that the critical line is remarkably linear in the 
(Li, L2) plane indicating that a simple sum Li + L2 controls the concentration of energy 
in a large part of parameter space. 



5 Spectrum of emitted gravitons 




Figure 1: Scattering of two homogeneous beams of size Lj = 1. Inclusive spectra 
|A;||/i(A;)p of emitted gravitons, as a function of \k\ close to criticality. Bottom to top: 
R = 0.44, 0.45, 0.46, 0.47, with = 0.470673. The two curves for each R are for n = 60 
and n = 70. 

Let us recall, from that the graviton spectrum is determined in terms of h{k) 
through: 



"^"^ GsR^\h{k)\\ (32) 



In FigH] we show a suitably normalized transverse-energy distribution of gravitons, 
\k\h{kY ~ as generated from the iterations (|26l) . The two adjacent curves for each R 
give an idea about the residual dependence on the "volume" n. As usual, the convergence 
with n is slower in the vicinity of the critical point, but the Figure suggests that the 
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behaviour near is rather regular. In particular there is no indication for a buildup of 
any divergence in the spectrum as R ^ R~ . 




0.5 1 1.5 2 2.5 3 



Figure 2: Same as in the previous figure but in configuration space. The various curves 
(brom bottom up) correspond to R = 0.45, 0.46, 0.47, 0.4706, 0.47064, 0.47065, and Rc = 
0.470673. 

A similar conclusion follows from Figl2] where the density profile of the gravitational 
field (7ri?)^r/i(r)^ = r(l — p{r)Y in transverse distance r is shown. The numerical cal- 
culations were done in configuration space following [Tj. There is no "finite volume", n, 
in this case and we can concentrate on the dependence on R. Again, no singularity in r 
develops as R ^ Rc and we can define a smooth limiting distribution at the critical point 
R~ . Note however, that the dependence on R close to Rc is rather strong (the change 
of R is tiny for the three uppermost curves), suggesting that the limiting distribution is 
attained with a large, possibly infinite, derivative. 

Figure [3] illustrates yet better this point. It turns out that the maximum value of the 
above density depends on {Rc — i?) as a simple square root. When left unconstrained, the 
best fitted power was always within ±1% of 1/2, remaining very stable against adding or 
removing initial/final data points. The solid line in Fig. [3] shows the fit where the power 
was actually fixed to 1/2. Similarly, we have found that the width of the distribution 
also behaves as Ci + C2{Rc — RY^"^, with finite coefficients Ci,C2. Therefore, indeed, the 
limiting distribution exists and is approached with an infinite derivative w.r.t. R. At 
the transition point, gravitons are emitted preferentially from half the distance from the 
source's edges. 

In [1] we have already analyzed the total multiplicity of emitted gravitons since it is 
related to the derivative of the action with respect to R. We found that it approaches a 
finite constant at R^ with a square root branch point at Rc (see figure 5 of p^j), with a 
"best fit" given by 0.138 — 0.46 (i?c — -R)°'^^^. This is hardly surprising given the above 
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results for the differential distribution. 




To close the circle we have also compared the Fourier transform of the p-space solution 
h{k) with h{r) obtained directly from the solution p(r) in x-space. The two agree locally 
within 2-3 %, the discrepancy being caused again by the finite discretization in momentum 
space. This is reassuring, not only in confirming again the consistency of the whole 
procedure, but also because, due to the infrared behaviour, transition between x- and p- 
representations is rather subtle. In particular, we find (c/. again Fig. [T])that the emission 
amplitude h{k) (and also Pi{k)) diverges at small momenta as 

hik) ~ i= , (33) 

and is exponentially damped at large k. Consequently, the action Q) is IR divergent and, 
even if the divergence is physically irrelevant (c/. the infinite Coulomb phase), has to be 
treated with care numerically. 

The large k behaviour of the spectrum can be qualitatively assessed from Fig. HI The 
distribution resembles much more an exponential than the gaussian shape of the sources 
we have put in (this is the C case of Sect. 3). Actually, a rather interesting structure 
emerges, as shown in that Figure. Let us use units in which L = 1 and increase R (i.e. 
the energy), starting from very small values. In that perturbative region the spectrum 
appears to consist of two exponentials separated by a "knee", i.e. the slope at small k, 
is smaller than the one at large k. However, as we increase R towards its critical value, 
the knee tends to disappear leaving behind an almost perfect exponential exp{—b\k\). 
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The slope b of the exponential (determined mainly by L in the perturbative region) now 
strongly depends on (R — Rc). Fig JS] illustrates the points we just made and indicates that 
the slope also has a rather singular behaviour, possibly of the form 6 ~ ci + C2(-R — -Rc)^- 








2 . 5 



5 



7 . 5 



10 



12 . 5 



15 



17 . 5 



k 



Figure 4: A logarithmic plot of the distribution ([32]) for = .10, .20, .40, 60, .80, .82, .83 
(bottom to top). 

One can finally try to determine the total multiplicity by integrating the spectra (either 
in a; or in p space). At small R/L one finds: 



Given that the average transverse energy is 0(1/ L) this corresponds to an average total 
transverse energy in the emitted gravitons given by: 



However, as one approaches Rc, this quantity clearly becomes 0{y/s) possibly implying 
that the transverse energy becomes a good estimate of the total radiated energy but 
also, unfortunately, that imposing energy-conservation, i.e. taking into account the back- 
reaction on the sources, becomes mandatory. 

6 Beyond the critical point: a first attempt 

There is obviously much interest in extending the numerical analysis into a BH phase. 
Although precise solutions for R > R^ are not available at this stage, we can get an 
approximate picture of their complicated structure using the momentum-space approach. 



(A^) ~ Gs{R/Lf . 



(34) 




(35) 
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R 



Figure 5: R dependence of the slopes from the previous figure (determined locally from 
10 < A; < 12.5). 

As already explained in Sect. 4, the discretized equations of motion fl22l) can be solved 
exactly (though numerically) as an algebraic system. This can be done for arbitrary i?, 
without any use of the iterative approach. The price is, of course, that we deal with 
a non-linear system of 3n equations for 3n variables, which is quite challenging even 
numerically. Restricting to identical sources reduces a problem to n variables. Still there 
are 3" solutions, and Mathematica has to generate all of them before we can choose 
the physically acceptable ones. To produce all 729 solutions for n = 6 takes about 2.5 
hours hence, in practice, the method is limited at present to n < 5. Still, even such 
a crude discretization reasonably reproduces the main features of, say, the momentum 
distribution. Therefore we go ahead with the simplified problem here and look for complex 
solutions above Rc leaving the, obviously possible, refinements for the future. 

Figure [6] shows the behaviour of two complex trajectorie^ as R moves from the disper- 
sive to the BH phase. We show only a discrete series of points with varying step AR (see 
below) to give also an idea of the speed (in R) along a trajectory. For n = 5, the critical 
point is at Rc = 0.3193. We begin by following the left trajectory (also reproduced by the 
recursion in this region) deep inside the dispersive phase (at R = 0.15): it moves towards 
its critical value Xc along the real expected. Checking the stability matrix at the 

fixed point we verified that, indeed, this is the only solution which is stable against the 
iterations fl2B]) . The second solution (which starts above Xc in the Figure) is also real, 
but unstable. We have chosen it to be the one that matches the recursive solution at 
Xc- The sampling points are separated by AR = 0.01 at the beginning of the recursive 
solution and then, as we approach Rc, the coarse graining was reduced to AR = 0.001. 
One readily sees that while approaching Xc the variation w.r.t. R increases in agreement 
with the findings of Sect. 5. For the unstable trajectory AR was set to 0.001 from the 

^By a trajectory we mean here a path traced by h{k), at one value of k, while changing R. 
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-0.1 
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Figure 6: Two complex trajectories h{k) (at k 
through the critical point. 



2.) parametrized by R moving 



beginning since we started tracing it much closer to Re- Again, its velocity increases as 
we go towards the critical point. As soon as one goes above Rc both trajectories acquire 
(complex conjugate) imaginary pats y{R)- Closer inspection shows that these imaginary 
parts grow like a square root, ^/x — Xc - typically for a threshold behaviour. Since both 
trajectories are complex conjugate to each other. Fig. [6] implies that the real part of 
the left trajectory continues to increase while that of the right trajectory changes from 
decreasing to increasing as we go through R^. Interestingly, the trajectories reveal a rich 
structure even when we go deeper inside the BH phase. We have followed them up to 
R = 0.65 changing the R-resolution: AR = 0.001 0.01 0.005 0.01 and finally to 
AR = 0.003 along the last segments. In these regions the trajectories appear to saturate. 
It remains to be seen if all these detailed features, like bending or saturation, are generic 
or are just artefacts of our small value of n. For example, we know that increasing n 
would move Rc by more than 30 %. 

However, one thing is clear: for R > Rc classical solutions of ( 126|) (and consequently 
also the on shell effective action) develop imaginary parts. This is turn implies a new 
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absorption of the elastic amplitude (on top of the one due to graviton production) calling 
for the opening of some new channel. The whole mechanism is somewhat reminiscent 
of the classic discussion of the "decay of the false vacuum" by Coleman and de Luccia 
[12] as a tunnelling process described through the contribution of complex saddle points 
to the functional integral. Use of similar ideas in this new context is presently under 
investigation [T3] . 

Let us conclude by stressing that the non-linearities captured by equations fl2B]) are 
essential for the above instabilities to occur. The gravitational attraction alone exists 
already in the lowest Born diagram, but it is not sufficient to produce the non-linearities 
of the metric that are essential for the buildup of CTS. Our results confirm that, instead, 
the class of diagrams selected in [2] appear to be sufficient for bringing out the emergence 
of such phenomena. 
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